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Abstract. The application of Fabry-Perot interferometers (FPIs) to the 
study of upper atmosphere thermodynamics has largely been restricted 
by the very low light levels in the terrestrial airglow as well as the limited 
range in wavelength of photomultiplier tube (PMT) technology Durinq 
the past decade, the development of the scientific grade charge-coupled 
device (CCD) has progressed to the stage in which this detector has 
become the logical replacement for the PMT. Small fast microcomputers 
have made it possible to “upgrade” our remote field sites with bare CCDs 
and not only retain the previous capabilities of the existing FPIs but ex- 
pand the data coverage in both temporal and wavelength domains. The 
problems encountered and the solutions applied to the deployment of a 
bare CCD, with data acquisition and image reduction techniques, are 
discussed. Sample geophysical data determined from the FPI fringe pro- 
files are shown for our stations at Peach Mountain, Michigan, and Wat- 
son Lake, Yukon Territory. 
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1 Introduction 

One goal within the National Science Foundation's Coupling, 
Energetics, and Dynamics of Atmospheric Regions Program 
(NSF-CEDAR) is to develop and deploy state-of-the-art in- 
strumentation at field sites scattered across the globe that will 
perform unattended automatic data acquisition of terrestrial 
airglow and auroral emissions. These data sets provide a 
means to characterize the weather at high-altitude regions of 
the atmosphere: the mesosphere, the lower and upper thermo- 
sphere, and the exosphere. Simple optical experiments con- 
trolled by small data storage instruments are one technique 
employed by aeronomers. In particular. Fabry-Perot inter- 
ferometers (FPIs), suitably outfitted, are capable of moni- 
toring neutral and ion temperatures, neutral horizontal and 
vertical winds, ion drifts, and composition measurements 
from stable environmentally controlled platforms. Other in- 
struments used by experimental optical aeronomers are dif- 
fraction grating spectrometers or spectrographs, spectropho- 
tometers, all-sky imaging cameras, and Fourier transform 
spectrometers. The common connecting thread of all these 
techniques is the desire to perform very low light level ex- 
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periments on temporally varying phenomena to extract geo- 
physical parameters with the lowest possible error bars. 

We describe the current technique used by a group of 
aeronomers at the University of Michigan to perform Fabry- 
Perot interferometry at remote field sites. These locations are 
typically difficult to reach, have little if any on-site support, 
are in harsh environments, and usually are isolated from the 
principal investigators for large parts of the observing season. 
The FPI - measures the line profile of an emission originating 
from a long-lived state of either atomic oxygen (ground or 
singly ionized state) or hydroxyl. The excited state has a 
lifetime sufficient for it to reach local thermodynamic equi- 
librium and, hence, provides a measure of the local temper- 
ature. Shifts in the peak location provide an indication of the 
line-of-sight velocity, which may easily be converted into 
either a horizontal (meridional and/or zonal) or a vertical 
wind component. The key measurement of these FPI airglow 
experiments is the 2-D fringe profile — intensity versus phase 
difference. We discuss the method employed to image the 
fringe profile and to extract the fringe profile from the image, 
as well as the difficult task of deploying a bare CCD as the 
sensor for an FPI airglow experiment. 

The FPI instruments deployed are all located in the 
northern hemisphere: Thule Air Base, Greenland; Sdndre 
Strdmfjord, Greenland; Watson Lake, Yukon Territory, Can- 
ada; and Peach Mountain, Michigan. The Greenland sites 
currently use an image plane detector 1 to gather light, while 
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the latter two employ a bare CCD as a detector. As we learn 
to deal with cryogen coolant, the S^ndre Strtfimtjord facility 
will be upgraded with a bare CCD. In addition, an all-sky 
imaging experiment was redeployed at Peach Mountain with 
a new upgraded bare CCD detector during the fall of 1992. 

2 Experiment 

The absolute radiative quantity measured, known as the in- 
tegrated intensity or surface brightness, is usually assumed 
to" be isotropic, and within the aeronomical community is 
scaled by the absolute unit, 4 the rayleigh (R), where I R is 
equivalent to Kf photons cm' 2 (column) s ' 1 . The absolute 
emission rates for the features observed by the FP1 are in the 
neighborhood of 10 to 100 R. in which 1 kR is typical of the 
limit of naked eye visibility. Hence, the FP1 airglow exper- 
iment must be capable of measuring very low light levels 
with minimal contribution from inherent noise sources. 

An example of our typical FPI airglow experiment is given 
in Meriwether et a). 5 All stations have basically the same 
physical layout, but use one of two differently styled detec- 
tors The image plane detectors employed in Greenland con- 
sist of only 12 channels, while the CCDs at Peach Mountain 
and in the Yukon have thousands of individual pixels that 
may be ordered into ~ 100 channels. Consequently, with ap- 
propriate optical elements, a CCD-based FPI may be designed 
to image several fringes simultaneously, an inherent multi- 
plex advantage over the Greenland image plane detector 
(IPD)-based FPI. 

3 Charge-Coupled Device Selection 

The primary requirements for the choice ot a bare CCD tot 
our FPI airglow experiments are low inherent thermal noise 
and low readout noise. The portion ot the CCD used in the 
FPI airglow experiment is a circular area closely aligned with 
the center of the detector and, hence, defects around the edges 
of the CCD are not important for our photometric work. This 
latter point permits some reduction in sensor cost, because a 
perfect detector is not required tor the experiment. Our ex- 
periments are conducted with visible- or near-infrared radia- 
tion, therefore, front-illuminated devices have satisfactory 
quantum efficiency for our requirements. In the last several 
years, our group has experimented with four different CCDs 
for use within FPI and all-sky imaging experiments: (l)a 
Fairchild CCD222 device, (2) a Thomson TH7882 CDA de- 
vice. (3) a Ford/Photometrics PM5 1 2 device, and (4) a Kodak 
KAF-14(X) device. The numerical order matches our chron- 
ological purchasing and. coincidentally, cost (the more recent 
the purchase, the greater the cost). Note that (2) and (3) are 
no longer in production. Device (1) was initially used with 
a lab bench FPI before installation in an all-sky imaging 
experiment. We have recently replaced the Fairchild chip 
with another Kodak device and have fielded this upgrade. 
Devices (2), (3). and (4) have been fielded in Ann Arbor, 
Peach Mountain, and Watson Lake, though we will soon 
move device (3) to Sdndre Strcimfjord. returning (2) to active 
use. Table 1 describes the devices according to manufacturer 
specifications. We have performed some checks on these 
specifications, though a report is usually issued on receipt ot 
each camera system. 

The Fairchild CCD is a fairly inexpensive design that has 
been used in several other aeronomical applications: for ex- 


Table 1 Specifications of bare CCD devices employed for FPI stud- 
ies. The primary characteristic that determines the usefulness of the 
CCD is the noise level. In use, the Thomson chip is binned 3x3, 
the PM512 chip is binned 4x4, and the Kodak chip is binned 8x8. 
This on-chip binning increases the dark noise component of the 
super pixel by the square root of the number of component pixels, 
but the readout noise associated with the digital number correspond- 
ing to the super pixel remains as listed. 


Device 

Number of 

Pixel 

Readout noise 

Dark noise 


pixels 

size (iim) 

(c- RMS) 

(e /sec pixel) 

Fairchild CCD222 

488 x 380 

12 x 18 

60 

18000 @ 24 °C 

Thomson TH7882 CDA 

384 x 576 

23x23 

1 1.5 w/ bias 

4.1 e/min @ -120 °C 

PM512 

516x516 

20 x 20 

8 

6.0 @ 45 °C 

Kodak KAF1400 

1317 x 1035 

6.8 x 6.8 

11.2 

0.003 @ -42 °C 


ample, see Torr et al. 6 The system was purchased in a single- 
detector CAM/CCD5000 design in which the CCD resided 
in the CAM5000 sense head with a Pel tier-effect cooler ca- 
pable of reducing the CCD die temperature by 2()°C below 
ambient. The control unit, a CCU5000 camera subsystem, 
provided output in RS170A television format, implying a 
1 ,4-V peak-to-peak 75-11 source impedance, composite an- 
alog video output. This translates into a 525-line. 30-Hz frame 
rate which is acceptable for bright signals but useless for 
low light level sources. Modifications to the hardware were 
made in our lab to (1) remove an IR filter from the sense 
head enclosure to extend the wavelength sensitivity ot the 
device, (2) permit time exposures with integrations ranging 
up to several minutes, and (3) further cool the sense head by 
placing it within a Products for Research (PfR ) thermoelectric 
cooler. The latter modification allowed the CCD to be cooled 
to about — 35°C, which reduced a very large thermal noise 
component. The final design was fielded at Thule Air Base 
within an all-sky imager devoted to monitoring auroral and 
airglow patches during the long polar winter nights. Data are 
stored in analog form on VHS videotape, though a frame 
crabber board within the data control system permits digi- 
tization of either real-time analog signals or previously re- 
corded images. 

Initial dark tests with the Fairchild CCD indicated that 
without the PfR cooler, the CCD was unacceptable for low 
light level studies. The detector saturated in 20 s, as indicated 
by the average 8-bit digital pixel value. Further test results 
with the CAM5000 placed inside the PfR cooler permitted 
observations of several minutes. However, airglow obser- 
vations were still difficult to make, though auroral measure- 
ments at Thule Air Base were possible, proving our concep- 
tual design. 7 

The Thomson CCD system was purchased from Photo- 
metries, Ltd., with a liquid nitrogen cooling requirement. The 
use of cryogen permits a drastic reduction in the operating 
die temperature of the CCD, with a consequent reduction in 
thermal noise. The system was ordered w ith a CC 200 camera 
electronics unit, a CC200 camera controller, and a side- 
looking CH2 10 cryogenic camera head. Note that these com- 
ponents, as described, are no longer available from Photo- 
metries, including the Thomson TH7882 CDA sensor. The 
entire experiment was initially controlled via a GPIB interface 
with a Micro Vax II controller, though later, an IBM Al 
clone supplanted the larger computer. Our first successful 
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airglow FPI observations were performed with this system 

fooi lab u m Ann Arb ° r ’ Michi § an - Beginning in September 
' , ' th,s system began remote unattended operation at a 
held site at a dark location on Peach Mountain. Michigan. 
Figure 1 contains two images taken with this system- 
fa) shows an image of an OI (5577-A) fringe, while (b) shows 
an image of the OH-Meinel (7,3) P,(3) line at 8920 A. 

The PM512 CCD system was also purchased from Pho- 
tometries, Ltd., with components similar to the Thomson 
system, though the CH210dewar was of the bottom-looking 
variant. This required the addition of another small mirror in 
the optical path to bend the light beam onto the sensor rather 
than placing the dewar on its side (LN2 would have dribbled 
out). The original set of software control routines designed 
to operate the Thomson-based system required some revision 
to handle the different pixel arrangement of the PM5P 
Otherwise, the PM512 and the Thomson systems were iden- 
tical. This system has been evaluated for long-term unat- 
tended operation at a field site on Peach Mountain, Michigan, 
and will soon be transferred to an aeronomical observatory 
in Sondre Strdmfjord, Greenland. Geophysical data were ac- 
quired in support of a CEDAR-sponsored mid-latitude, me- 
sosphere. and lower thermosphere thermodynamic project, 
as well as ground truth support for the orbiting (JARS 
mission. 

The Watson Lake, Yukon, facility was deployed in No- 
vember 1991 as a remote upper thermosphere monitoring 
stanon. The FPI was an old design originally deployed in 
Alaska and then moved to Calgary, Alberta, in 1980 The 
detector used in Calgary was an IPD device that acquired 
lower thermosphere measurements 5 during the early 1980s. 
Equipment failure and a lack of funding precluded any further 


experimental measurements until 1991. The optical hardware 
was returned to Michigan and a bare CCD-based system was 
purchased from Photometries, Ltd. This time, a Kodak 
KAFI400 CCD was chosen as the sensor. The rest of the 
system comprised a CH250 camera head, an LC200 liquid 
circulation unit, and an IBM AT interface card. This system 
differed markedly from the older Photometries versions in 
that it was liquid cooled rather than cryogen cooled, thus 
making the detector head more compact. The Kodak CCD 
operates solely in the multipinned phase (MPP) mode 8 and, 
thus, the expected enhanced thermal noise contribution from 
operating at -40°C is much reduced. Figures 1(c) and 1(d) 

^. ay ‘"" g ! S acquircd b - v this s y s,cm: I (c) shows an image 
of OI (6300-A) fringes, while 1(d) shows an image of HeNe 
laser fringes. 

The chronological history of our attempts at using and 
upgrading very low light level hardware shows that we^have 
been very successful. Three FPIs use a bare CCD device and 
an all-sky imager has been outfitted with another Kodak de- 
vice. The sensitivity range of the CCD is such that near-IR 
radiation may be easily observed, adding a new dimension 
to our experiments — hydroxyl emissions become brighter as 
one moves from the visible to the near IR, allowing ground- 
based remote observation of mesospheric altitudes. As funds 
become available, remaining FPIs and diffraction grating 
spectrometers will be upgraded with the versatile CCD. 

4 Fabry-Perot Interferometry with a Bare 
Charge-Coupled Device- 

4.1 Fitting a Circular Fringe onto a Square Matrix 

Initial observations of frequency-stabilized HeNe laser fringe 


OPTICAL ENGINEERING / February 1994 / Vol. 33 No. 2 / 459 







Fie 2 Initial observations of the transmitted circular fringe pattern from a HeNe-illuminated Fabry* 
pilot interferometer (a) the original uncorrected fringe pattern, (b) a dark count integration (c) a flat 
S!ma e of a whi.eligh, source, and (d) the suitably ** M **** 

and aspect ratio corrections for the rectangular pixel size of the Fairchild CCD. 


patterns were performed with the Fairchild CCD system. 
Figure 2 displays one of the first observations. Flere a white- 
frosted globe illuminated by the HeNe laser was used to 
scatter the laser radiation to evenly illuminate a 2-in.-diam 
set of etalon plates. The interference pattern was then focused 
onto the Fairchild CCD by means of a simple lens. The image 
was digitized with an 8-bit frame grabber board, stored, and 
then displayed on a monitor attached to a Micro Vax II sys- 
tem. Figure 2(a) shows an original raw image, while a sub- 
sequent dark image is shown in Fig. 2(b). A flat-field image 
of a tungsten bulb illuminating the globe is shown in Fig. 
2(c). These last two images are required for our image pro- 
cessing: the dark frame is removed from both the FleNe and 
the flat-field images, and then the dark-corrected HeNe image 


is divided by the dark-corrected flat field. The first technical 
problem that we encountered was caused by the rectangular 
size of the pixels within the Fairchild CCD. In addition, the 
CCD222 was developed with an interline transfer architecture 
in which adjacent vertical photosites are covered by opaque 
aluminum shields. After each integration, a charge is trans- 
ferred from the exposed vertical columns to the adjacent 
shielded columns so that another integration may take place 
while the previous one is read. The opaque shield is 
18X18 pm. while the active photosites are 12X 18 p.m. 
Hence, the initial images will never show a circular pattern 
when one is input, as observed. However, by applying an 
aspect ratio correction, a circular image may be formed as 
shown in Fig. 2(c). The utility of the Fairchild device, though, 
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is limited by the opaque shields, because one-half of the fringe 
pattern will never be observed. This and the fact that the 
CCD222 quickly saturated with thermal noise made it un- 
acceptable tor airglow FPI observations. 

The next set of test patterns was acquired with the 
Thomson-based CCD system. HeNe laser fringe profiles were 
of adequate quality and attempts were made to reduce radially 
the fringe pattern from a 3-D nature (two linear and one 
intensity dimensions) to the conventional intensity versus 
phase difference scheme. To measure geophysical parameters 
accurately from the resultant fringe profile requires a pro- 
cedure to perform this reduction. 9 We know that an ideal 
transmitted pattern from a FPI is represented by an Airy 
function that has intensity maxima at integral values of the 
order of interference in, given by 

in\ = 2\it cosf) , 

where \ is the wavelength of the incident radiation, t is the 
gap between the etalon plates. 0 is the angle of the light path 
through the etalon relative to the normal, and p. is the re- 
fractive index of the material or gas between the plates. The 
input angular dependence may be replaced by the relationship 
between radial distance from the center of the fringe pattern 
r and the effective local length of the system/,, 

0 = tan \r/fj , 

and because the angles are small, may be expanded in a Taylor 
series yielding 



Thus, the wavelength is linear with respect to the square of 
the linear distance from the center of a transmitted FPI fringe 
pattern. 

Figure 3 displays a result from a simulation performed to 
test the sensitivity of the r 2 reduction technique to an accurate 
description for the centerof the FPI pattern. Figure 3(a) shows 
the initial Airy profile as a function of radial distance using 
the appropriate instrumental values for the Thomson-based 
CCD system. The fringe profiles clearly bunch together as 
the phase difference term passes through several orders of 
interference. Here, the area covered by each physical pixel 
was further subdivided into 100 equal area subpixels. The 
effect of integrating the intensity levels of subminiature pixels 
onto actual size pixels is shown by Fig. 3(b). The hatched 
line indicates the original input Airy function, while the his- 
togram display indicates what the CCD would actually ob- 
serve. Applying the software aperture technique of Mulli- 
gan, 9 a simulated r 2 version was generated and is shown in 
Fig. 3(c). A slight broadening occurs in the data reduction 
technique resulting in a reduction of peak intensity for dif- 
ferent orders. 

The next step in the simulation was to offset the center 
of the fringe pattern and to redo the r reduction. Figure 4(a) 
displays the simulated fringe pattern. Various runs were per 
formed with a sample result shown in Fig. 4(b). Here the 
value of the center coordinate for one axis has been offset 
by one-half of a pixel dimension, resulting in a fringe profile 
indicated by the solid line. The dashed line shows the ex- 





Fig. 3 (a) Airy profile simulation of a transmitted fringe pattern on 
illumination by a frequency-stabilized HeNe laser. Appropriate val- 
ues for reflectivity, gap spacing, and the illuminating wavelength 
were included in the numerical simulation. The figure shows the ex- 
pected input to a bare CCD. (b) The histogram function shows the 
simulated intensity pattern across a central axis of a bare CCD on 
illumination. The hatched line shows the input function, (c) The result 
of applying the r 2 reduction technique described in the text. The 
locations of maximum intensity are now evenly spaced as a function 
of bin number. 
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(b) 


Fig. 4 (a) The input fringe pattern and (b) the result of mismatching 
the location of the center of (a). Here, the center was misplaced by 
one-half of 1 pixel, resulting in a broadened and doubled i 2 profile, 
compared to the anticipated profile shown as a dashed line. 


peeted profile for perfect alignment. Not only are the higher 
order fringes unnecessarily broadened, but a double peak 
appears because the simple r reduction program cannot han- 
dle elliptical patterns. Typical HeNe finesse for the central 
fringe is of the order of 15, but with poor alignment, outer 
fringes show finesse values of 5 or less. Consequently, a great 
deal of energy is expended during start-up operations to en- 
sure proper alignment of the optical elements for the best 
possible setup prior to remote unattended operation. 

4,2 Bias and Noise Reduction 

A difficulty encountered with low light level airglow obser- 
vations is degradation caused by charge transfer efficiency 


(CTE) problems, that is, losses due to inefficient transfer of 
charge from one photosite to the next as a CCD is read. 
Among various others, 10 the Thomson TH7882 CDA sensor 
was known to exhibit a reduction in transfer efficiency when 
the number of electrons per pixel was small (—50 e ). The 
technique employed by Photometries to impiove the CTE 
was to illuminate the CCD weakly and uniformly with a 
preflash, or optical fat zero, equivalent to 100 e /pixel, prior 
to any observation. As shown in Janesick et al., 1 the CTE 
is dramatically improved when this technique is used. How- 
ever, inherent in the addition of this bias charge to the electron 
count per pixel is an increase in the noise floor of the device 
due to the shot noise that the additional illumination provides. 
This noise is in addition to the other usual sources of noise: 
shot noise due to illumination by airglow or aurora, dark 
current generation due to charges of thermal origin, and 
preamplifier or readout noise. The normal techniques em- 
ployed to improve the final image are to address the reduction 
of these noise sources by increasing the integration period to 
reduce the shot noise, to reduce the die temperature of the 
sensor to eliminate dark current generation, and to bin pixels 
on chip so that readout noise is only added to the * ‘super 
pixel" rather than to the individual pixels. 

Our observations of oxygen and hydroxyl airglow at mid- 
latitudes, samples displayed in Fig. 1, are usually performed 
with a 5-min integration period. Initially, observations were 
performed such that the system gain was determined by the 
ratio of the full well capacity of each super pixel and the full 
range of the analog-to-digital converter. It was soon discov- 
ered that unnecessary digitization noise was being added to 
the final images — the difference between maximum intensity 
and minimum intensity values tor airglow fringe profiles was 
of the order of unit values for typical profiles, in which the 
minimum value included the sky plus all noise components. 
Experiments were performed to optimize system gain tor the 
greatest expected airglow signal values. In practical terms, 
this meant reducing the value for the analog-to-digital sat- 
uration limit to a much smaller number. This implied that 
the full well capacity was greatly underutilized and that shot 
noise components contained higher digital values. Histo- 
grams of 5-min dark count (plus bias) observations indicated 
that when a representative bias frame was removed from the 
dark frame, the average digital value was zero (ignoring ran- 
dom cosmic ray strikes and other “hot-pixel values). Hence, 
dark count observations were no longer performed during 
sky observations to increase the observing frequency of air- 
glow data, though dark frames may still be required for longer 
integrations. 

The final sky frame still includes a contribution from the 
bias preflash. As shown in Janesick et al., l() CTE performance 
in a CCD responds dramatically to die temperature. At the 
optimum system gain, the contribution of the bias component 
to the sky frame was still substantial and required a process 
to carefully estimate and remove it from the frame. A high 
SNR bias frame is usually acquired by averaging 81 separate 
bias frames prior to performing sky observations. The cry- 
ogenic camera heads that we use are temperature stabilized 
to within it 0. 1°C, but at the system gains that are employed, 
we have found that average bias frame values slowly drift 
during long-term nighttime observations. This behavior is 
associated with small temperature drifts immediately follow- 
ing the delivery of cryogen to the detector dewar, amounting 
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to roughly 5% of the average nighttime bias value (note that 
the CCD is not operated at the cryogen temperature, but 
usually at - IOO°C). The routine used to remove the bias 
component from the sky frame utilizes this high SNR bias 
frame by subtracting this amount from the sky frame. We 
have found this to be adequate except during intervals when 
the dewar housing the camera head is being refilled with 
liquid nitrogen. Our operational answer to this problem was 
to control the refill timing to coincide with dusk and dawn, 
the latter to ensure that the detector dewar would remain cool 
during daylight periods. We note here that with the new 
Kodak KAF 1400 based system in Watson Lake, preflashing 
of the sensor is not used because airglow emissions are in- 
herently brighter (due to the high latitude) and the etalon 
plates have a much larger area (6-in. diameter versus 4-in. 
diameter) resulting in higher count rates. The Watson Lake 
airglow sample shown in Fig. 1(c) was acquired in 1 min. 

4.3 Geophysical Observations 

Two stations are currently operating with bare CCD sensors: 
Watson Lake, Yukon, using the Kodak chip, and Peach 
Mountain, Michigan, using the Thomson device, though the 
PM512 chip has been thoroughly tested at the latter site. 
Figure 5 displays a sample from the first stage of processing 
the raw fringe data acquired in Watson Lake on day 60, 1992, 
Here the FPI was observing emission from the OI (6300- A) 
feature in five different directions: vertical and at a 45-deg 
angle in the four cardinal directions. Both frames show the 
1 reduction results for the five measurements per cycle. The 
integration period per direction is determined by first ac- 
quiring a sample 10-s image. Watson Lake is sufficiently 

close to the northern auroral oval that the possibility of auroral 

contamination is quite high, and as such could saturate the 
image, hence, a quick integration is necessary per direction 
from which a time scaling factor is determined for the final 
recorded image. The range available is from the original 10 
to 1 80s, determined primarily from a numerical simulation 
for anticipated uncertainty estimates for the line-of-sight 
winds and temperatures that may be extracted from the mea- 
surement, here, about 1 0 m/s and 25 K, respectively. The two 
sets of five fringe profiles clearly show the wavelength shifts 

that the bare CCD FPI combi nation i s capable of measuring 

here amounting to about 200 m/s in the meridional direction. 

The Peach Mountain station routinely observes emissions 
from the Meinel band of hydroxyl. This feature provides a 
means to measure the neutral winds in the mesosphere. Fringe 
profiles similar to those shown in Fig. 5 have been routinely 
acquired for the last several years. A data reduction technique 
has been developed derived from that described in Killeen 
and Hays. 1 1 Instrumental calibration is performed by pressure 
scanning the interferometer while observing a frequency- 
stabilized HeNe laser. A set of 160 images is acquired pro- 
viding a scan over two full free spectral ranges for each of 
the r annular bins (in our analysis, we chose 150 such bins). 
The periodic behavior of the fringes is numerically reduced 
to Fourier coefficients, which are shifted in wavelength to 
the appropriate sky wavelength, here 8920A. Sky observa- 
tions also contain instrumental contributions convolved with 
the actual Doppler-broadened skyline profile. The fitting pro- 
cedure chosen is to assume a Gaussian profile for the airglow 
emission that contains four unknowns: ( I ) width of the profile 
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Fig. 5 Radially reduced OI (6300-A) airglow fringes acquired at Wat- 
son Lake, Yukon. Five fringes are shown in each panel. Shifts in 
peak location are evidence of strong horizontal neutral wind com- 
ponents. 


that is proportional to the kinetic temperature of the emitting 
species, (2) wavelength shift of the center of gravity of the 
profile that is proportional to the line-of-sight Doppler ve- 
locity, (3) a scaling factor that is proportional to the absolute 
brightness of the emission feature, and (4) a remnant back- 
ground that is proportional to weak sky continuum. This 
analytical representation of the airglow emission line profile 
is then convolved with the numerical representation of the 
instrumental profile at the appropriate wavelength and then 
compared with the actual observed r~ profile. Better guesses, 
in a least squares sense, are then made for the four unknowns 
until convergence is reached. 12 Figure 6 displays horizontal 
neutral wind results for two successive periods 13 in 1990. 
The meridional (N, S) and the zonal (E, W) components are 
shown with error bars as a function of local time, with N and 
E being positive. Both low- and high-frequency temporal 
variations are evident in the wind component data, and these 
are current topics of study within the aeronomical commu- 
nity. 

5 Conclusions 

The current set of scientific grade CCDs has been shown to 
be capable of performing aeronomical quality observations 
of airglow emissions with Fabry-Perot interferometers. Both 
routine (OI 6300-A) and unique (Meinel OH 7,3 band) ob- 
servations have been acquired with several different bare 
CCDs. The fringe quality is satisfactory and standard reduc- 
tion techniques are able to extract high-quality geophysical 
data from the observations. The CCD is a versatile sensor 


OPTICAL ENGINEERING / February 1994 / Vol. 33 No. 2 / 463 





NICIEJEWSKI, KILLEEN, and TURNBULL 


Day number 259/260. 1990 — RAB CCD FP1 



Day number 260/261, 1990 RAB CCD FPI 



22 24 2 4 6 8 10 

Local Time (hr) 


22 24 2 4 6 

Local Time (hr) 


Fin fi Wind determination from fringe profiles of the Meinel (7,3) band of hydroxyl as observed at a 
mid latitude^ s^te^T^o^^ccessive nights of horizontal neutral wind data are displayed as meridional 
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and, with proper precautions, may be used in an unmanned 
automatic data-taking mode with little or no maintenance. 

Several initial problems have been addressed and solutions 
have been found: (1) the use of full-frame transfer devices 
with square pixels is far superior to the interline transfer 
devices, (2) cooling of the sensor and the use of the MPP 
mode are satisfactory for the integrations required for the 
brighter airglow emission features, (3) careful attention to 
the choice for the center of a circular fringe pattern is nec- 
essary to reduce unnecessary broadening of the r fringe 
profile and to keep the useful finesse as high as possible, 
(4) chips that require an optical fat zero or bias contribution 
must be carefullv characterized for low light level observa- 
tions, and (5) standard data reduction techniques may be eas- 
ily applied to data acquired by bare CCD-based FPls. 

The bare CCD has quickly revolutionized the method of 
performing unattended remote field site observations of the 
thermodynamics of atmospheric regions. The detector is ro- 
bust and has an excellent quantum efficiency over the spectral 
regions of interest, a satisfactory dynamic range, and a lin- 
earity in response at the count rates acquired. Our group is 
planning to complete the upgrade of the set of field sites that 
we maintain and is currently considering further novel mod- 
ifications to lab-based CCD FPI systems, such as imaging 
fringes from two different frequencies simultaneously on the 
same device and imaging fringes onto concentric ring CCDs. 
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